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Depending on the rule for tree growth, the forest-fire model shows either self-organized criticality 
with rule-dependent exponents, or synchronization, or an intermediate behavior. This is shown 
analytically for the one-dimensional system, but holds evidently also in higher dimensions. 
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During the past years, systems which exhibit self- 
organized criticality (SOC) have attracted much atten- 
tion, since they might explain part of the abundance of 
fractal structures in nature Q. Their common features 
are slow driving or energy input (e.g. dropping of sand 
grains, increase of strain, growing of trees) and rare dis- 
sipation events which are instantaneous on the time scale 
of driving (e.g. avalanches Q, earthquakes 0, fires 
In the stationary state, the size distribution of dissipation 
events obeys a power law, irrespective of initial conditions 
and without the need to fine-tune parameters. 

While the critical behavior of the sandpile model is rel- 
atively well understood, it is in general not clear, under 
which conditions SOC behavior occurs. Models of the 
above-mentioned type show a variety of different phe- 
nomena, depending on details of the model rules. When 
e.g. the form of the driving force in the earthquake model 
|| is modified, a transition from SOC to synchronization 
(activity bursts covering the whole system) is found || . 
When in the same model the boundary conditions are 
changed from open to periodic, the SOC behavior is re- 
placed by a periodic state j5|] . Another earthquake model, 
which includes inertia, shows a superposition of a power 
law and a peak at large avalanches, where most of the 
energy is released Q. A continuous forest-fire model |t]] 
shows periodic behavior or finite avalanches, depending 
on the value of a parameter || . 

Another question in the context of SOC deals with the 
universality of the critical exponents. The critical be- 
havior of the sandpile model seems to be robust with 
respect to various changes, and the critical exponents of 
the forest-fire model do not depend on the lattice symme- 
try ||1J or the presence of a further parameter . In 
contrast, the critical exponents of the earthquake model 
depend continuously on the degree of energy conserva- 
tion A similar dependence of critical exponents 
on the model parameter is found in a new version of the 
forest-fire model with tree conservation fl3|| . 

Since a great part of the above-mentioned work is per- 
formed numerically, there is substantial need for analyt- 
ically understanding these various phenomena. In this 
paper, I study analytically a generalized version of the 
forest-fire model in one dimension. Different rules for tree 
growth lead to either SOC behavior with rule-dependent 



exponents, or to synchronization, or to an intermediate 
state with a superposition of a power law and a peak 
at large fires. These results are confirmed by computer 
simulations. Similar behavior is expected in higher di- 
mensions. 

The model is defined as follows: Each site in a one- 
dimensional system of length L is either occupied by a 
tree, or it is empty. Trees on neighboring sites belong 
to the same tree cluster. During one time step of size 
dt, each tree is struck by lightning with probability /. 
Trees struck by lightning and all other trees in the same 
cluster burn down and turn to empty sites immediately. 
Empty sites have a distribution of life times P(r), i.e. 
P(r)dt is the probability that a site which just burned 
down will remain empty exactly until the time t. The 
probability that a site which has been empty for a time 
t becomes occupied by a tree during the next time unit 
dt is therefore P(t)dt/ f. P(r)dr. It is reasonable to re- 
quire that the mean life time f of empty sites is finite: 
f = L P(r)TdT < oo. This means that P(t) decays 
faster than r~ 2 for large r. The condition that burning 
is fast compared to tree growth places the model in the 
class of systems with slow driving and instantaneous dis- 
sipation events. The mean number of trees growing per 
unit time is L(l — p)/f, and the mean number of light- 
ning strokes per unit time is Lfp, where p is the mean 
tree density in the stationary state. The mean number 
of trees destroyed per lightning stroke consequently is 

s = (l-p)/ffp. (1) 

When / becomes very small, the mean number of trees 
destroyed by a fire is large. If there are only large fires 
but no small ones, large parts of the system burn down 
simultaneously and may therefore be synchronized. If 
there are fires of all sizes up to a cutoff size, the system 
is close to a critical point and shows scaling over many 
orders of magnitude. In this paper, we always assume 
that the system size L is so large that no finite-size ef- 
fects occur. The cutoff is then a function of / and P(t) 
(see Eq.(||) below). 

In the original forest-fire model, each empty site be- 
comes occupied by a tree with probability p, and conse- 
quently P(t) = p exp(— pr). In the limit f/p — > 0, the 
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size distribution of forest clusters in the stationary state 
is then essentially a power law with a cutoff at sizes of the 
order (p//)/ln(p//) Q. In the following, we will deter- 
mine analytically the size distribution of forest clusters 
for different P(r). 

To this purpose, let us consider a string of k neigh- 
boring sites in the system. If the lightning probability / 
is sufficiently small, lightning does not strike this string 
before all its trees are grown. Starting with a completely 
empty state, the string passes therefore through a cycle. 
Trees grow on the string, until it is completely occupied 
by trees. Then the forest in the neighborhood of the 
string will also be quite dense. The forest on the string 
is part of a forest cluster which is much larger than k. 
Eventually that cluster becomes so large that it is struck 
by lightning with a nonvanishing probability. Then the 
forest cluster burns down, and the string again becomes 
completely empty, and the cycle restarts. The mean time 
T(fc), which it takes for k trees to grow on k empty sites, 
satisfies 

/>oo 

/ P(r)dT cx l/k. (2) 
JT{k) 

This equation can e.g. be obtained from the condition 
that the tree density is larger than 1 — l/k for T = T(k). 
A more rigorous derivation gives the same result. We 
assumed that the probability that lightning strikes the 
string during that time is negligeable. The maximum 
string size for which our considerations are valid is there- 
fore given by the condition 

/•I'M 

fcmax/ / (T(fc max ) - T) P{T)dT < 1. (3) 

Jo 

After time t, an initially empty string of k < /c max sites 
is occupied by m(t) = k J* P(r)dr trees on an average. 
For not too small values of m, fluctuations around this 
mean value are relatively small. The probability Vkijn) 
to find m trees on a string of k < fc max sites therefore 
satisfies for 1 << m < k 

V k (m) cc [{dm/dt^J- 1 = (kPit^))- 1 , (4) 

where t m k is the time after which the initially empty 
string is occupied by m trees. Let n(s) be the number 
of clusters of s trees, divided by the number of sites L. 
n(s) is identical to the probability that a string of s + 2 
neighboring sites is occupied by s neighboring trees, with 
one empty site at each end, 

n(s) =P S+2 ( S )/( S + 2 ) ^2P s+2 ( s )/s 2 . (5) 

For small s, there are deviations from this law due to 
the discreteness of the lattice. Remember that the result 
Eq.(||) is valid only for s < s max = fc max — 2. For larger 
s, clusters are struck by lightning with a nonvanishing 
probability before they can grow to larger clusters, and 
therefore Eq.(||) does not hold any more. 



The mean number of trees destroyed per lightning 
stroke can be expressed in terms of the cluster size dis- 
tribution, 

OO Smax OO 

s = ^2s 2 n(s)/p= s 2 n(s)/p+ ^ s 2 n(s)/p. (6) 

1 1 s max + l 

Here, we have to distinguish two different situations: De- 
pending on the precise form of n(s), the second term can 
be neglected in the limit / — » with respect to the first 
term, or it dominates. In the first case, the system has 
fires of any size and therefore is in a critical state, in the 
second case, the dynamics are dominated by large fires. 
The first case will occur if P(t) has a long-time tail. Then 
the neighborhood of the tree which is struck by lightning 
is likely to contain surviving empty sites which eventu- 
ally stop a fire. The forest density approaches the value 
1 in the limit / — ► 0, since otherwise s could not diverge. 
The second case will occur if tree growth is finished after 
a finite time. Then large regions of the system will be- 
come simultaneously occupied by a dense forest and burn 
down together. In this case the mean forest density will 
be smaller than 1. 

To make these statements more precise, let us specify 
the different possible cases: 

(i) P(t) has a long-time tail, i.e. P(t) oc t~ q for large 
r, with a > 2. Sites which have been empty for a time 
t become occupied by a tree during the next time step 
with probability (a — V)dt/t, where a is now a parame- 
ter and not an exponent. From Eqs.(||) and (||) we find 
TO) cx fcV(o-l) and s max oc /-!+!/« With Eqs.(g) and 
(||) we obtain the size distribution of forest clusters 

n(s) oc a -2+V("-i). 

The forest density p = Y^Li sn i s ) cannot be larger than 
1, and consequently n(s) assumes a normalization factor 




FIG. 1. Scaled size distribution sn(s)// 1 of fires as 
function of s/s max - The parameters are L — 10000 and a = 3, 
for 1// = 100000 (solid), 50000 (dotted), 25000 (dashed), 
12500 (long dashed), 6250 (dot-dashed), 3125 (solid), 1562.5 
(dotted). The straight line has the slope —0.5. 
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With Eq.(0) we find s oc s max , which leads with Eq.(|l]) 
to 1 — p oc /V". These results represent SOC behavior, 
and the values of the exponents depend on a, i.e. they 
are non universal. The analytic results are confirmed by 
simulations. Fig.[l] shows the collapsed curves for the size 
distribution sn(s) of fires for L = 10000 and a = 3, for 
different values of/, confirming the relations for n(s), for 
Smax, and the normalization factor oc f x ^ a . Fig.|| shows 
the density of empty sites as function of /, confirming 
the relation 1 — p oc J 1 /". 
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FIG. 2. Density of empty sites as function of 1/f for 
L — 10000 and a — 3. The solid line has the slope —1/3. 

(ii) P(r) decays exponentially fast: P{t) ~ exp(— t@) 
for large r. This situation comprises also the original 
forest-fire model. The exponents are obtained from case 
(i) by taking the limit a — > oc. In addition, there oc- 
cur logarithmic corrections. With Eqs.(0), (^), and 

we find T(k) ~ (In(k)) 1 ^, w(Mw)) 1//3 oc 1/f, and 
n(s) oc s~ 2 /ln(l//). For (3=1, these results have al- 
ready been analytically derived and confirmed by simu- 
lations in [p) . 

(iii) All trees grow within a time To, i.e. P(t > To) = 
and J T ° P(t)cIt = 1. In this case /c max is essentially given 
by fcmax fx 1//T . The large-scale dynamics of this sys- 
tem can be obtained by performing a scale transforma- 
tion: fc max sites form together a big site which becomes 
occupied exactly after To steps, and the lightning prob- 
ability for this coarse-grained system is fc max / oc 1/Tq. 
This is a model with deterministic tree growth, which 
shows synchronization and will be considered in the next 
paragraph. Before, let us take a look at the short- 
scale dynamics: The size distribution of clusters smaller 
than s max is obtained in the same way as before. If 
lim T ^T P(t) / 0, oc, we find n(s) oc s~ 3 and conse- 
quently Y^* max s 2 n(s)/p oc ln(l//To). Since we know 
(Eq.(P) that s oc 1//T , we conclude from Eq.(|) that 
most fires are large and are not described by the power 
law sn(s) oc s~ 2 . Fig.^ shows the size distribution of fires 
for P(t) = 1/T , T = 100, L = 10000, and / = 0.00001. 
In addition to the mentioned power law, there is a peak 



at large fires, indicating that large segments of the system 
burn down together. 
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FIG. 3. Normalized size distribution of fires for 
P(t) = 1/100, T = 100, L = 10000, / = 0.00001. The 
smooth line is the theoretical result Eq.(5), which is valid for 

The special case lim T ^o P(Tq — r) oc r Q for a > or 
-1< a < leads to n(s) oc s -3+a/(i+«). 

(iv) Last, let us consider the case of deterministic tree 
growth P(t) = S(Tq — t), which is also the coarse-grained 
version of situation (iii). Here it is easy to understand 
how synchronization arises: If two neighboring sites hap- 
pen to be both occupied by a tree, they burn down during 
the same fire. Consequently, both trees regrow simul- 
taneously and burn down simultaneously for all future 
times. So the system consists essentially of synchronized 
blocks which turn to trees and burn down simultaneously. 
Neighboring blocks join, if the block which grows first is 
not struck by lightning before its neighbor grows. So the 
number of blocks decreases with time, and their mean 
size increases. In the stationary state, the mean block 
size is so large that each block is struck by lightning 
immediately after it grows, so that neighboring blocks 
cannot join any more. This stationary state is periodic 
with a period To. In the limit /To — > 0, the block size 
becomes very large, and eventually the whole system is 
synchronized and fires every To time steps. Fig.^ shows 
the number of blocks as function of time for L = 10000, 
/ = 0.005, and To = 50, starting with a random state. 

Generalization to higher dimensions is straightforward: 
It is clear that the case P(t) = 5(Tq—t) leads to synchro- 
nization of large blocks of the system in any dimension. 
On the other hand, it is known that the system is SOC 
for P(t) = p exp(-pt). Therefore, it can be expected 
that when varying P(t) the critical exponents change 
their values, and that finally the system becomes syn- 
chronized. However, the calculations performed in this 
paper cannot be applied to higher dimensions. In con- 
trast one dimension, fire can now burn regions which are 
not completely dense, and part of the trees in that re- 
gion survive the fire, thus generating correlations which 
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are not present in the one-dimensional model. 

To summarize, I have shown that depending on the 
rule for tree growth the stationary state of the forest-fire 
model shows either SOC or synchronization, or a super- 
position of both. In the SOC state, the critical expo- 
nents vary continuously when P(t) is changed. These re- 
sults are similar to those for the earthquake model which 
shows SOC or synchronization depending on the form of 
the driving force Qj. In it is shown that the SOC 
behavior in the earthquake model is related to the ten- 
dency of neighboring sites to synchronize. Neighboring 
sites which topple during the same avalanche are likely 
to topple also the next time during the same avalanche. 
The forest-fire model shows a similar feature in the SOC 
state: A region which has just been crossed by a fire is 
empty (in one dimension) or almost empty (in two di- 
mensions). This region is with a high probability again 
covered by an almost dense forest in the moment when 
it is ignited. So neighboring sites which burned down 
during the same fire are likely to burn down simultane- 
ously also the next time. Due to the randomness of tree 
growth, this "synchronization" is not perfect. The degree 
of synchronization increases when tree growth becomes 
more deterministic, first showing a peak at large fires, 
and finally perfect synchronization of large clusters. 



The results of this paper might have implications for 
excitable media (e.g. spreading of diseases, autocat- 
alytic chemical reactions, propagation of electrical activ- 
ity in neurons or heart muscles). These systems essen- 
tially have three states which are called quiescent, ex- 
cited, and refractory. If there is an excitation, it spreads 
to the quiescent neighbors. After excitation, a region is 
refractory to further stimulation and needs some time to 
recover its quiescent state. The model discussed in this 
paper describes an excitable medium in the limit where 
excitation spreads fast on the scale of the refractory time, 
and where spontaneous excitation is rare. In cases where 
this limit can be realized experimentally, such a system 
should show either synchronization or SOC, depending 
on whether the transition from the refractory to the qui- 
escent state and from the quiescent to the excited state 
are deterministic or stochastic. 

I thank S. Clar and one of the referees for their com- 
ments on the manuscript. This work was supported by 
the Deutsche Forschungsgemeinschaft (DFG) under Con- 
tract No. Dr 300/1-1, and by the NSF grant No. DMR- 
93-03667. 
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FIG. 4. Number of synchronized blocks as function of time 
for L = 10000, / = 0.005, and To = 50. 

When the stochastic rule for lightning strokes is re- 
placed by a more deterministic one, the model also be- 
comes more synchronized. Lightning might strike only 
trees which have a certain age, or only clusters which 
have a certain size. In these cases, there are only large 
fires in the stationary state and no small ones. However, 
the size distribution of forest clusters smaller than a cut- 
off size is not considerably affected by these changes since 
Eq. (|j) was derived under the only assumption that light- 
ning does not strike a string of size k < fc max , before all 
its trees are grown. 
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